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CENTRAL LIMIT THEOREM UNDER VARIANCE UNCERTAINTY 


DMITRY B. ROKHLIN 


Abstract. We prove the central limit theorem (CLT) for a sequence of independent 
zero-mean random variables £j, perturbed by predictable multiplicative factors A j 
with values in intervals [ A,. Ay] . It is assumed that the sequences A j, A j are bounded 
and satisfy some stabilization condition. Under the classical Lindeberg condition 
we show that the CLT limit, corresponding to a “worst” sequence A j, is described 
by the solution v of one-dimensional G-heat equation. The main part of the proof 
follows Peng’s approach to the CLT under sublinear expectations, and utilizes Holder 
regularity properties of v. Under the lack of such properties, we use the technique of 
half-relaxed limits from the theory of viscosity solutions. 


1. Introduction 


Consider a sequence of independent one-dimensional random variables with 

zero means and finite variances crj = > 0. Put = ]U™ =1 cr|, £ > 0 and assume 

that the Lindeberg condition 

1 n 

Ln{z) — — X! E (^/{|g 7 i>es»}) 0; n y oo (1.1) 

s n j=1 


is satisfied. Then, by the classical central limit theorem (CLT), for any bounded 
continuous function / : 1 1 we have 


lint E/ 

n—>oo 



E/K), 


( 1 . 2 ) 


where ( has the standard normal law. 

In this paper we assume that the variance of is not known exactly and may belong 
to an interval. Our goal is to obtain the “least upper bound” for the quantity (II. 2ft 
under such model uncertainty. The result, as well as its proof, are similar to those 
obtained by Peng mm and the followers (TTJ [2D1 H] under the nonlinear expectations 
theory paradigm. It appears that can be described in terms of the solution v of 
a nonlinear parabolic equation, called G-heat equation. One of the objectives of the 
present paper is to show that this description also comes from a classical problem 
statement, and need not be linked to the nonlinear expectations theory. 
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To give a precise problem formulation, consider a filtered probability space 

(^,p, mr= o) 

and an adapted sequence of random variables such that Efj = 0, E £j = of G 

(0, oo) and ^ is independent from Let (Aj)°L 0 be an adapted sequence, whose 

elements Aj take values in deterministic intervals [A -, A^], 0 < A ? < A j. Considering 
the sequence rjj = one can regard the multipliers Aj_i as a “predictable per¬ 

turbation” of the original sequence The intervals [A^cr,-, \j_\aj\ indicate possible 
standard deviations of rjj. 

Assumption 1. The Lindeberg condition (11.11) is satisfied. 


Assumption 2. The sequence A j is bounded by a constant A. 


Assumption 3. The sequences A j; A j satisfy the following stabilization condition-. 

n —1 2 

= + ^°, n ^°° ( L3 ) 
j=o 

for some A > A > 0. 


Put Bj = [A 2 , A 2 j\, B = [A 2 , A 2 ] and denote by 

d H (Bj, B ) = max{|A 2 - A 2 1, |A 2 - A 2 |} 

the Hausdorff distance between these intervals (see, e.g., [H, Chapter 2]). Condition 
(11.31) is equivalent to the following one: 


n ~ l a 2 

Y, -y-d H (Bj, B ) ->• 0, n -» oo. 

3 =0 Sn 


(1.4) 


In the summability theory the transformation 

+ _ Pl a l + ‘ ‘ + Priori 
tn — : ; 

Pi H-E Pn 


Pn> 0 


of a sequence (cq)^ is called a Riesz mean (see [HU Section 1.4], [5J Section 3.2]). 
By Assumption [3l the sequence dn(Bj,B) is summable to 0 by the Riesz method, 
determined by the sequence pi = of. Furthermore, the Lindeberg condition implies the 
Feller condition 


Oj 

lim max — = 0 

n—too 1 <j<n S n 


(1.5) 


(see, e.g., [3], Chapter 6, §28). In particular, s n —> oo. Hence, the Riesz summation 
method, defined above, is regular (see [ 161, Theorem 1.4.4]), and if dn(Bj, B ) —> 0, then 
the Assumption [3] is satisfied. We also mention a necessary and sufficient condition 
for (11.41) to hold true, given in [4] (Lemma 3.2.14). This result is applicable since 
s 2 n 0 by m- 
Note that from the identity 


n —1 


£ 



1 
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it easily follows that 


lim 

n—> oo 


n —1 

X 

3=0 


a: 


j+ it 2 _ t 2 
c 2 ^ ’ 


lim 


n—1 

X 

3=0 


a 


3+ 1 \ 2 _ ->2 

„2 —j — ' 


Denote by 21 n the set of adapted sequences Ag = (A j)] =0 with values in [A^, Aj]. Our 
goal is to describe the quantity 

AA+iJ , (1-6) 

which can be loosely characterized as the least upper bound of (II. 2p under variance 
uncertainty. 

The main role in this description is played by the solution of the nonlinear parabolic 
equation 

v t + \ sup_ (. \ 2 v xx ) =v t + \ (a\A - A 2 ^) = 0, (t, x ) G [0,1) x M, (1.7) 

z ag[a,a] 2 v 7 

satisfying the boundary condition 

v(l,x)=f(x), igR. (1.8) 

In the context of the CLT under sublinear expectations, equation (I1.7K appeared in 
ra- It was called G-heat equation in ra- As is mentioned in [5], such equation arises 
in various applications in control theory, mechanics, combustion, biology, and finance. 
It is known also as a Barenblatt equation: see, e.g., [9]. 

One can obtain (JTTTj) by considering A j as a control sequence, writing down dynamic 
programming equations for discrete time finite horizon optimization problems 



and passing to the limit as n — » oo. This approach was proposed in [T8] in the case of 
identically distributed (multidimensional) random variables £j. However, in the present 
context, it seems that this method requires hypotheses, which are stronger than the 
Lindeberg condition. Thus, we follow Peng’s approach, which takes equation (H77j) 
as a starting point, and utilizes a deep result on the existence of its solution in an 
appropriate Holder class. 

Put Q = [0,1] x M, 

||fr||0;R = SUp|/l(x)|, ||0||o;Q= Slip \g(t,x) |, 

ccgm (t,x)eQ 


f 1 "- 1 

= lim sup E / — 

MOO AJ- ka-1 V Sn f~0 


[9\a;Q 


[^]q;R = Slip 


Xly^X2 


\h(xi) - h(x 2 )| 


xi - x 2 


a ’ 


a e (0,1], 


sup 

(t-i j 


\g(ti,xi) -g(t 2 ,x 2 )| 

(In - nl 1/2 + \xi - x 2 \) a 


a e (0,1], 
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and consider the Holder spaces C 2+ "(IR), C l+a ^ 2,2+a {Q) with the norms 

||^||c 2 +“(R) = ||^||o;K + || hx || 0;R + ||h ra ||o ; R + [h xx ] q; r, 

IM|c , 1 +“/ 2 ' 2 +“(Q) — IM|o ; Q + ||<7®||o ; Q + ||<7t||();Q + ||<7®®||o;Q + [, 9t\a-,Q + [dxx]a-,Q- 
Under the assumptions / G C ,2+ “(M), a G (0,1]; A > 0 the existence of a classical 
solution v G C 1+a '/ 2,2+a '(Q) (with some of a' G (0,1]) of (II.7ft . (JUS]) was proved by- 
Krylov: see [TO] (Theorem 1.1 or Theorem 5.3). 

If A = 0 then only the existence of a viscosity solution is guaranteed. Let us recall 
this result along with related definitions. Put Q° = [0,1) xl and assume that / is 
a bounded continuous function: / G Ct(M). A bounded upper semi continuous (use) 
function v : Q H > R is called a viscosity subsolution of (11.71) , (11. 8p if 

v(l,x) < f(x), igR, (1.9) 

and for any (t,x) G Q° and any test function p G C 2 (M 2 ) such that (t,x) G Q° is a 
strict local maximum point of v — p on Q°, the inequality 

- p t {t,x) - ]- sup_ (A 2 p xx (t,x)) < 0 (1.10) 

1 Ae[A,A] 


holds true. To define a viscosity supersolution, one should consider a bounded lower 
semicontinuous (lsc) function v, a strict local minimum point of v — p, and reverse the 
inequalities (II.9ft . (11.101) . 

We will use the following comparison result. Consider a viscosity subsolution u and 
a viscosity supersolution w of (1 1.70 , (11.81) . Since we require (11.71) to be satished in the 
viscosity sense at the lower boundary of Q, by the accessibility theorem of |5], we have 


u(0,x)= lirnsup u(t,y ); 

(t,y)G(0,l)xR, 
t^0,y—>x 


w(0,x) = liminf wit,u) 
(t,y)e(0,l)xR, 
t—>0,y^x 


and by the comparison result of [7] (Theorem 1) it follows that u < w on Q. 

A bounded continuous function v : Q i —> M is called a viscosity solution of (11.71) . 
(11.81) . if it is viscosity sub- and supersolution. The existence of a continuous viscosity 
solution of (11.71) . (j 1.8 p for / G C),(K) is well known from the theory of optimal control. 
The stochastic control representation of such solution can be found in [T9] (Chap. 4, 
Theorem 5.2). 


Theorem 1. Let f be a bounded continuous function, and let v be the continuous 
viscosity solution of \1. 1\) , \ 1.8\) . Then, under Assumptions^^, we have 2z? = n(0,0). 

It is interesting to compare Theorem [T] with related results obtained in the framework 
of sublinear expectations theory. Besides the original result of Peng mm, which is 
discussed in [TS], we mention the papers HUM, where the random variables were 
not assumed to be identically distributed. We will discuss only the result of [2D], which 
extends HU- The result of JSj concerns the multidimensional case. 

Let us briefly describe the construction of a sublinear expectation space (9,77,E), 
which allows to rewrite the expression (11.61) in terms of a sublinear expectation. This 
construction is, in fact, the same as in [18, Section 4], where some more details are 
given. Consider the space of sequences 9 = {{yi)°L 1 : k G R}, and introduce the space 
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of random variables % as follows: % = U^'Hn, where V. n is some linear space (we do 
not go into details) of functions Y = i/j(yi,... ,y n ) of n variables. Define the sublinear 
expectation by the formula 

EY = sup EV>(A 0 6M> ■ ■ •, A„_i£ n /cx n ). (1.11) 

A o _le 2 t o _1 

Let Yi be the projection mappings: Y+y ) = y t . One can show that Y n is independent 
from (Ij,..., E n _i) in the sense of sublinear expectations theory (see (15], Definition 
3.10). By (II.lip we get the following representation for J?f: 


2^ 


lim sup E / 

rwoo A”- 1 e2i"-i 



limE flT-Yj 


d =1 


Let us apply Theorem 3.1 of [2D] to the sequence Y t . We have 
I(±Af = sup _ E(±Aj_i^/c7j) = 0, 

— — 1 di — l] 

Ilf = sup _ E(Aj_i£j/cu) 2 = X-i, 

Ai_ iG[A i _i,Ai—i] 

-E(-lf) = - sup _ E (-(Ah 6M) 2 ) = At,- 

A«—iG[A i _i,Aj_i] 


Besides a condition, identical to Assumption [3] in [20j it is assumed that 

I|lf 2+5 = sup _ E = A^E|^j/crj| 2+5 < M, (1.12) 

Ai—i G [A i —i ,Aj_i] 


lim V 

n .—Von ' ^ 


3 =1 


CD 


2+5 


= 0 


(1.13) 


for some M > 0. S > 0. Note that (I1.12p was used in PH in this form, although it was 
not clearly formulated (see condition (3) of Theorem 3.1 in [2D]). The result of [20j 
tells us that YY = E/(Z), where Z is a G'-normal random variable with 

G( s ) = 2 ( s+ ^ — s A 2 )- 

By the characterization of the G-normal distribution (see, e.g., [15j, Example 1.13) this 
is equivalent to the assertion of Theorem [0 

Thus, under the assumptions (II . 12[) . (II . 13[) (instead of Assumptions [T] [2]) , Theorem 
[1] follows from the result of |20j . It is easy to see that (I1.12p implies Assumption [2} 

< fE(A^|6/^| 2+5 )' V(2+5) 


-2 .2/ 2^ 1/2 


Aj-i = (E(K-iZtdt) 
and (11. 12ft . (II. 13ft imply that 


1 


AEDA-ky 


3 = 1 


l {\Zj\>£Sn}) < £ 5 s 2+8 


E 5 

3 =1 




2+5 


M 


o2+<5 


E 

3 = 1 


2+5 


o, 


n —> oo. 
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The last condition is slightly weaker than Assumption [U and coincides with the latter 
if lim inf j^oo A; > 0 . 

Note that if there is no model uncertainty: A ? = A j = 1, then Theorem |T] reduces to 
the classical CLT, mentioned at the beginning of the present paper. This is not the case 
with the result of [ 20 j, since in this case the conditions ( 1 1 . 12 j) . (1 1.13 [1 are stronger then 
the Lindeberg condition. We also mention that |20j deals only with classical solutions 
of the G-heat equation, so the case A = 0 is, in fact, not considered there. However, the 
sublinear expectations theory is able to handle the degenerate case via perturbation 
methods, see [Hj (the proof of Theorem 5.1), [B] (the proof of Theorem 3.1). 


2. Proof of Theorem [H 

(i) We first consider the case / G G 2 +Q (M), a > 0 and A > 0. Put 

A 


J ol 


Xj + \ — Xj -f— -£j+ i, j — 0,..., n — 1, A"o — 0 ; tj — -£■ 

Sn k =o Sn 

Since the solution v of (11.7k (ll.Sji belongs to v G cT+“7 2 . 2 +«' (Q) , we can apply Taylor’s 
formula: 

n—l 

u(l, X n ) - w(0,0) = ( V ( t J+U X 3 + 1 ) - v{tj,X j+ 1 ) + v{tj, x j+ 1 ) - v{tj, x j)) 

3=0 

n - 1 / 1 \ 

= ^ ^ ( v t(tj, X j+i){tj+i ~ tj) + v x (tj, Xj)(Xj _|_i — Xj) + —v xx (tj,Xj)(Xj -|_i — Xj) j , 

3=0 A ' 

where tj = tj+/3(tj + i—tj), Xj = Xj + r y(Xj + i~Xj), /3 ,7 G [0,1]. By the independence 
of Xj and Q + \ we conclude that E(v x (tj,Xj)(Xj +1 — Xj)) = 0. Thus, 


n—l ^2 

Eu(l,X n )-u(0,0) = E^ 


3=0 


a „ f ^ A 2 £ 2 - 

v t{tj,Xj + 1 ) + -X- J 2 + v xx (tj , Xj) ) = J n + I n , 
1 a j +1 


n—l 2 


Jn ~ E ^ 


cr: 


l+i 


3=0 


v t{tj, Xj) + 


\ 2 c2 
A ! Sl +1 
9 2 ‘'zaU 

z a i+i 


n-l 2 


fe.Y) =eE 


cr; 


l+i 


l=o 


A 2 , 

Vt T ~^~v xx ) (tj, Xj), 


In = 


n—l 0 
07 


i=o 


E E -X v *(5v A Vi) - **(*!> *l) + 


2 ^ + i 


Ai (tj , Aj) v xx (tj, Xj )) 


We can rewrite J n as + J 2 , where 


n—l 2 


• 1 » = e E 


a: 


l+i 


l=o 


A + ^ ( A V* - A 2 Wx ) ) (tj , Afj), 


1 ^a 2 

t 2 _ _ — ^ ^ 

^ n r\ 

3=0 ~ n 


n-i 9 

e E A 1 (A - Wi + (a 2 - a>-) (%, Aj), 
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From the definition of v we see that J\ = 0. Furthermore, from the stabilization 
condition (II .3j) it follows that 


l_^a 2 


1+1 <(\) - A 2 )«i + (A 2 - ) fe, X s ) < CM n 0, 


sup + - p e V 
z i=0 

^ oo, since the second derivative of v is uniformly bounded. On the other hand, 
choosing a sequence 

Aj = ^jl{v x ^xqx)} - Aj^ M (q,Y 3 )<0} ; j > 1, 
with an arbitrary Aq, we get an opposite inequality 


1 n ~ l rr 2 
sup J 2 n >- 
Aj’hr - 1 Z j =0 b n 


K - A )v+ + (A - A )u- fo, x,) > —CM n -A 0. 


Combining all these results, we conclude that 


lim sup J n — 0. 


n )-oo « ji —i 


(2.14) 


AJ _1 ea n - 1 


Now consider I n — I x n + 7 2 + 7^: 

n-! a 2 


! n = E E dF W+ Y+i) - +C A',)), 


5=0 


Ai ^ ^ 2 o ( V xx{tjjXj) V xx (tj,Xj) \ I{ |£j + 1 |>es„}, 

,y=o 2 

£2 y2 ^ 

In = ^ ^ o ~o“ ( v xx{tj,Xj) — V xx (tj, Xj)j -f{|£ i+ i|<es n }- 
j=o Z 

By the Holder continuity of Vt we have 


I III 


71—1 

<cte£ 

5=0 


a: 


5+1 


I A? 


la '/ 2 


+ |X 


5+1 


Ah- 



Using the inequality E|<^^ / • + l| Q, < (E£j +1 ) Q,//2 = and the independence of \ 3 and 

^j+i, we obtain the estimate 



From (II.5|1 it follows that 7* —> 0. 

Furthermore, since the sequence A j is bounded and the second derivative of v is 
uniformly bounded, by the Lindeberg condition we get 

|7 2 | + CL n {£) —y 0, n —y oo. 
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The last term 1% is estimated with the use of the Holder continuity property of v xx : 


rc-l c2 \ 2 

lhl<cE^h±iJ- 

3=0 n 

Therefore, 


A? 10 +i I 


- 1 ^2 


'{hj+l|<£Sn} — 


< CA 2+a ' - J i L£ °' = CA 2+a 'e a '. 


3=0 Sn 


lim sup \I n \ =0. 


(2.15) 


A£ _ 1 e 2 D 


From (12.141) and (12.15p it follows that 

2 z? = lim sup Ef(X n ) = lim sup Eu(l, X n ) — u(0, 0). 


Ag €21" 


A ™” 1 € 21 " 


So, we have proved the theorem in the case / G C 2+a , A > 0. 
(ii) Now assume that A = 0. Put 


71—1 


X £ n = — J^(A 2 + e 2 ) 1/2 £ j+1 , Af £ =lim sup E f (X*). 

n^An-i^-x 

The intervals [/i., Jl^ = [(Aj+e 2 ) 1 / 2 , (A^ + e 2 ) 1 / 2 ] stabilize to [e, (A 2 + £ 2 ) 1 / 2 ] in the sense 
of Assumption [3j 

n— 1 9 

(l ^ 2 - + £2 )l + \£ ~ £ 2 l) ^ °> n ->■ °°- 

3=0 Sn 

By part (i) of the proof, we infer that Af £ = u e (0, 0), where v £ satisfies 

1 ( -2 r 2\f a ,e 5 + _ J2( e \-\ _ 


h £ + 2 (( r + £ )^) +_e ( t '») J =0 ’ xeQ°-, u £ (l ,x) = f(x), xeR (2.16) 

in the classical sense. Let v be the continuous viscosity solution of the limiting problem 


\ _ 2 

v t + 2^ v xx = °j xeQ°; v(l,x) = f(x), x G 
The desired result is a consequence of the relations 


(2.17) 


A? lim sup Ef(X n ) — limAf £ , v(0, 0) = lim u £ (0,0), (2.18) 


"^“Arha "- 1 


£—>0 


£—^0 


which we are going to prove. 

Since we still assume that / G C' 2+Q (M), this function is uniformly Lipschitz contin¬ 
uous. Put ip £ ( A) = (A 2 + e 2 ) 1 / 2 — A. We have 

c ( /»-i ' 2N 1/2 

|e !(K) - e/(a'„)| < ce|a; - x n \ < - e 


,3=0 


n-i a2 


1/2 




3=0 Sn 
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since sup A>0 ^ e (A) = e. Thus, 

2zf £ — Ce < liminf sup E f(X n ) < limsup sup E f(X n ) < J? £ + Ce, (2.19) 

n-xx> A n-l Ga „_]. n^oo 

limsup < liminf sup E f(X n ) < limsup sup E f(X n ) < liminf J/f e . 

£—>■0 ri— >0 ° AQ _1 e2t n_1 n—>oo A «-ig a 7i-i £—>0 

These estimates imply the first equality in (12.18j) . 

Furthermore, dehne the half-relaxed (or weak ) limits of v e by 

v(t,x)— liminf v £ (s,y), v(t, x) = limsup v £ (s,y), (t,x)(zQ. 

(s,y)—>(t,x), 

e-vO 

The function v (resp., y) is use (resp., lsc): see [2] (Chap. 5, Lemma 1.5). 

Take p G (7 2 (IR 2 ) and assume that z = (t, x) G Q is a strict local maximum point of 
v — p on Q. Then there exist sequences —* 0, z*, = ( t k ,Xk ) G Q such that Zk —$■ z, 
v £k (zk) —$■ v(z), and Zk is a local maximum point of v £k — <p on Q: see [2] (Chap. 5, 
Lemma 1.6). 

If t G [0,1), then G [0,1) for sufficiently large k and 

-<Pt{zk) ~ sup (A 2 <p xx {z k )) < 0, 

AG[efe,A+£fc] 

since v £k is a viscosity solution of (12.161) . Passing to the limit as £& —> 0, we get the 
inequality 

- <p t (z) - sup_(A 2 p xx (z)) < 0, (2.20) 

Ag[0,A] 

which means that v is a viscosity subsolution of (12 .1 7p on Q°. 

Let t — 1. If there are infinitely many tk < 1, then we again obtain (12.20[) as above. 
Moreover, we can change the test function p to ip = <^ + c(l —£), c > 0 since (1,3;) is still 
a strict local maximum point of v — ip. Substituting (p in (12.201) . we get a contradiction: 

c - <p t (z) - sup_ (A 2 (p xx (z)) < 0, for any c > 0. 

Ae[0,A] 

Thus, for sufficiently large k, we have v £k (zk) = f(xk) and v(z) = lim fc _ > , 00 f(xk) = f(x). 

We have proved that v is a viscosity subsolution of (I2.17|) . Similarly, one can prove 
that y is a viscosity supersolution of (12.171) . By the comparison result of [7], mentioned 
in Section [H we have v < y on Q. The converse inequality v > y is immediate from 
the definition. We infer that v = v = y is a continuous viscosity solution of (12.171) . and 
the second equality in 02.1 8p holds true: 

u(0, 0) < lim inf u £ (0, 0) < lim sup if (0, 0) < u(0, 0). 

This completes the proof of Theorem |T] in the case A = 0. 

(iii) It remains to consider the case / G C&(R). It is not difficult to show that there 
exists a function f £ G (^“(M) such that \ f(x) — f £ (x)\ < e: see, e.g., [IT] . Furthermore, 
consider a function x £ C°°, 

X(x) = 1, |a:| < 1; x(x) = 0, |x| > 2, 




















10 


DMITRY B. ROKHLIN 


and put g £ {x ) = x( £l ^ x )f e ( x )- We have 

|E f(X n ) - Eg £ (X n )\ < |E f(X n ) - Ef £ (X n )\ + |E f £ (X n ) - Eg £ {X n )\ 

n ~ l A 2 a 2 

<e + CP(e 1/2 \X n \ > 1) < £ + CeEXl < e + - J 4 ±i = (1 + Ck 2 )e. 

3=0 Sn 

From this estimate we obtain the inequalities of the form (12.1 9ft with 

2zf £ = lim sup E g £ (X n ). 

Just mentioned inequalities imply that 

JSf:=lim sup Ef(X n ) — limJf £ . (2.21) 

n ^°° e ->° 

Denote by V £ , the viscosity solution of (II.7J) , ('ll .8(1 . corresponding to the terminal 
condition g £ instead of /. Since g £ G C 2+ “(M), we have 

£> e = V E (0,0) (2.22) 

by the result, already proved. 

Finally, note, that the convergence g £ (x) = x( £l 2x )f £ ( x ) —*■ f ( x ), £ —* 0 is uniform 
on compact sets. It follows that 

iimiiif g e (y) = lim sup g £ (y) = f(x). 

y—y—tx 
£— s —^0 

Using this fact, by the method of half-relaxed limits, applied above, it is easy to prove 
that 

lim D £ (0,0) = u(0,0). (2.23) 

£—>0 

From (I2.2ip - (I2.23I) we conclude that JX = u(0,0). The proof of Theorem |T] is com¬ 
plete. 
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